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Abstract We investigate the constrained Minimal Super- 
symmetric Standard Model (cMSSM) in the light of con¬ 
straining experimental and observational data from pre¬ 
cision measurements, astrophysics, direct supersymmetry 
searches at the LHC and measurements of the properties 
of the Higgs boson, by means of a global fit using the pro¬ 
gram Fitting. As in previous studies, we find rather poor 
agreement of the best fit point with the global data. We also 
investigate the stability of the electro-weak vacuum in the 
preferred region of parameter space around the best fit point. 
We find that the vacuum is metastable, with a lifetime signif¬ 
icantly longer than the age of the Universe. For the first time 
in a global fit of supersymmetry, we employ a consistent 
methodology to evaluate the goodness-of-fit of the cMSSM 
in a frequentist approach by deriving p-values from large 
sets of toy experiments. We analyse analytically and quan¬ 
titatively the impact of the choice of the observable set on 
the p-value, and in particular its dilution when confronting 
the model with a large number of barely constraining mea¬ 
surements. Finally, for the preferred sets of observables, we 
obtain p-values for the cMSSM below 10%, i.e. we exclude 
the cMSSM as a model at the 90% confidence level. 
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1 Introduction 


Supersymmetric theories QH) offer a unique extension 
of the external symmetries of the Standard Model (SM) 
with spinorial generators Q. Due to the experimental con¬ 
straints on the supersymmetric masses, supersymmetry must 
be broken. Supersymmetry allows for the unification of the 
electromagnetic, weak and strong gauge couplings 6 §- 
Through radiative symmetry breaking 1^1^, it allows for 
a dynamical connection between supersymmetry breaking 
and the breaking of SU(2)xU(l), and thus a connection be¬ 
tween the unification scale and the electroweak scale. Fur¬ 
thermore, supersymmetry provides a solution to the fine- 
tuning problem of the SM ||^[T^, if at least some of the su¬ 
persymmetric particles have masses below or near the TeV 
scale GD- Furthermore, in supersymmetric models with 
/^-parity conservation 112 131, the lightest supersymmetric 
particle (LSP) is a promising candidate for the dark matter 
in the universe |fT4l[T5|. 


Of all the implementations of supersymmetry, there is 
one which has stood out throughout in phenomenological 
and experimental studies: The constrained Minimal Super- 
symmetric Standard Model (cMSSM) |T§[T7). As we will 
show in this paper, albeit it is a simple model with a great 
set of benefits over the SM, it has come under severe ex¬ 
perimental pressure. To explain and - for the first time - to 
quantify this pressure is the aim of this paper. 


The earliest phenomenological work on supersymme¬ 
try was performed almost 40 years ago p^[T3][T8]j20) in 
the framework of global supersymmetry. Due to the mass 
sum rule GD, realistic models require local supersymme¬ 
try, or supergravity 116 22 -^. The cMSSM is an effective 
parametrisation motivated by realistic supergravity models. 
Since we wish to critically investigate the viability of the 












2 


cMSSM in detail here, it is maybe in order to briefly recount 
some of its history. 

The cMSSM as we know it was first formulated in 1^ . 
However, it is based on a longer development in the con¬ 
struction of realistic supergravity models. A globally super- 
symmetric model with explicit soft supersymmetry break¬ 
ing added by hand was first introduced in p7) . It is 
formulated as an SU(5) gauge theory, but is otherwise al¬ 
ready very similar to the cMSSM, as we study it at collid¬ 
ers. It was however not motivated by a fundamental super¬ 
gravity theory. A first attempt at a realistic model of sponta¬ 
neously breaking local supersymmetry and communicating 
it with gravity mediation is given in At tree-level, it 
included only the soft breaking gaugino masses. The soft 
scalar masses were generated radiatively. The soft breaking 
masses for the scalars were first included in | [29l[^ . Here 
both the gauge symmetry and supersymmetry are broken 
spontaneously p4| . In pO) the first locally supersymmetric 
grand unified model was constructed. Connecting the break¬ 
ing of SU(2) X U( 1) to supersymmetry breaking was first pre¬ 
sented in Q, this included for the first time the hi- and tri- 
linear soft-breaking B and A terms. Radiative electroweak 
symmetry breaking was given in Q. A systematic presenta¬ 
tion of the low-energy effects of the spontaneous breaking 
of local supersymmetry, which is communicated to the ob¬ 
servable sector via gravity mediation is given in |31 32\ . 

Thus all the ingredients of the cMSSM, the five pa¬ 
rameters tanjS, sgn(/i), Aq were present and un¬ 

derstood in early 1982. Here Mq and M 1/2 are the com¬ 
mon scalar and gaugino masses, respectively, and Aq is a 
common trilinear coupling, all defined at the grand unified 
scale. The ratio of the two Higgs vacuum expectation val¬ 
ues is denoted by tanjS, and /r is the superpotential Higgs 
mass parameter. Depending on the model of supersymme¬ 
try breaking there were various relations between these pa¬ 
rameters. By the time of p5) , no obvious simple model of 
supersymmetry breaking had been found, and it was more 
appropriate to parametrise the possibilities for phenomeno¬ 
logical studies, in terms of these five parameters. In many 
papers the minimal supergravity model (mSUGRA) is often 
deemed synonymous with the cMSSM. However, more pre¬ 
cisely mSUGRA contains an additional relation between Aq 
and Mq reducing the number of parameters . 

The cMSSM is a very well-motivated, realistic and con¬ 
cise supersymmetric extension of the SM. Despite the small 
number of parameters, it can incorporate a wide range of 
phenomena. To find or to exclude this model has been the 
major quest for the majority of the experimental and phe¬ 
nomenological community working on supersymmetry over 
the last 25 years. 

In a series of FITTING analyses |34-^ we have con¬ 
fronted the cMSSM to precision observables, including in 
particular the anomalous magnetic moment of the muon. 


(g — 2)n, astrophysical observations like the direct dark mat¬ 
ter detection bounds and the dark matter relic density, and 
collider constraints, in particular from the LHC experiments, 
including the searches for supersymmetric particles and the 
mass of the Higgs boson. 

Amongst the previous work on understanding the 
cMSSM in terms of global analyses, there are both those 
applying frequentist statistics | |38] - |58) and Bayesian statis¬ 
tics | |59H70) . While the exact positions of the minima de¬ 
pend on the statistical interpretation, they agree on the over¬ 
all scale of the preferred parameter region. 

We found that the cMSSM does not provide a good de¬ 
scription of all observables. In particular, our best fit pre¬ 
dicted supersymmetric particle masses in the TeV range or 
above, i.e. possibly beyond the reach of current and future 
LHC searches. The precision observables like (g —2)^ or 
the branching ratio of B meson decay into muons, BR(Bj — 
fifi), were predicted very close to their SM value, and no 
signal for dark matter in direct and indirect searches was ex¬ 
pected in experiments conducted at present or in the near 
future. 

According to our analyses, the Higgs sector in the 
cMSSM consists of a light scalar Higgs boson with SM- 
like properties, and heavy scalar, pseudoscalar and charged 
Higgs bosons beyond the reach of current and future LHC 
searches. We also found that the LHC limits on supersym¬ 
metry and the large value of the light scalar Higgs mass 
drives the cMSSM into a region of parameter space with 
large fine tuning. See also 171 ■ ^ on fine-tuning. We thus 
concluded that the cMSSM has become rather unattrac¬ 
tive and dull, providing a bad description of experimen¬ 
tal observables like (g — 2)^ and predicting grim prospects 
for a discovery of supersymmetric particles in the near fu¬ 
ture | [76) . 

While our conclusions so far were based on a poor agree¬ 
ment of the best fit points with data, as expressed in a rather 
high ratio of the global to the number of degrees of free¬ 
dom, there has been no successful quantitative evaluation of 
the ’’poor agreement” in terms of a confidence level. Thus, 
the cMSSM could not be excluded in terms of frequentist 
statistics due to the lack of appropriate methods or the nu¬ 
merical applicability. 

Traditionally, a hypothesis test between two alternative 
hypotheses, based on a likelihood ratio, would be employed 
for such a task. An example for this is e.g. the search for the 
Higgs boson, where the SM without a kinematically acces¬ 
sible Higgs as a “null hypothesis” is compared to an alterna¬ 
tive hypothesis of a SM with a given accessible Higgs Boson 
mass. However, in the case employed here, there is a signif¬ 
icant problem with this approach: The SM does not have a 
dark matter candidate and thus is highly penalised by the 
observed cold dark matter content in the universe. (It is ac¬ 
tually excluded.) Thus, the likelihood ratio test will always 
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prefer the supersymmetric model with dark matter against 
the SM, no matter how bad the actual goodness-of-fit might 
be. 

Thus, in the absence of a viable null hypothesis without 
supersymmetry, in this paper we address this question by 
calculating the p-value from repeated fits to randomly gen¬ 
erated pseudo-measurements. The idea to do this has existed 
before (see e.g. izz))’ but due to the very high demand in 
CPU power, specific techniques for the re-interpretation of 
the parameter scan had to be developed to make such a result 
feasible for the first time. In addition to the previously em¬ 
ployed observables, here we included the measured Higgs 
boson signal strengths in detail. We find that the observed 
p-value depends sensitively on the precise choice of the set 
of observables. 

The calculation of a p-value allows us to quantitatively 
address the question, whether a non-trivial cMSSM can be 
distinguished from a cMSSM which, due to the decoupling 
nature of SUSY, effectively resembles the SM plus generic 
dark matter. 

The paper is organised as follows. In Sec.|^we describe 
the method of determining the p-value from pseudo mea¬ 
surements. The set of experimental observables included in 
the fit is presented in Sec.|^ The results of various fits with 
different sets of Higgs observables are discussed in Sec. 
Amongst the results presented here are also predictions for 
direct detection experiments of dark matter, and a first study 
of the vacuum stability of the cMSSM in the full area pre¬ 
ferred by the global fit. We conclude in Sec. 


2 Methods 

In this section, we describe the statistical methods employed 
in the fit. These include the scan of the parameter space, 
as well as the determination of the p-value. Both are non¬ 
trivial, because of the need for theoretically valid 

scan points in the cMSSM parameter space, where each 
point uses about 10 to 20 seconds of CPU time. Therefore, 
in this paper optimised scanning techniques are used, and 
a technique to re-interpret existing scans in pseudo experi¬ 
ments (or “toy studies”) is developed specifically for the task 
of determining the frequentist p-value of a SUSY model for 
the first time. 


2.1 Performing and Interpreting the Scan of the Parameter 
Space 

In this section, the specific Markov Chain Monte Carlo 
(MCMC) method used in the scan, the figure-of-merit used 
for the sampling, and the (re-)interpretation of the cMSSM 
parameter points in the scan is explained. 


2.1.1 Markov Chain Monte Carlo method 


The parameter space is sampled using a MCMC method 
based on the Metropolis-Hastings algorithm |78-^. At 
each tested point in the parameter space the model pre¬ 
dictions for all observables are calculated and compared to 
the measurements. The level of agreement between predic¬ 
tions and measurements is quantified by means of a total 
which in this case correspo nds to the “Large Set” of observ¬ 
ables introduced in Section Is.lIH 

X ~ (Omeas Opred) COV meas Opred) + Zlimits) (1) 

where Omeas is the vector of measurements, Opred the cor¬ 
responding vector of predictions, cov the covariance matrix 
including theoretical uncertainties and X]?mits sum of all 
X^ contributions from the employed limits, i.e. the quantities 
for which bounds, but no measurements are applied. 

After the calculation of the total at the n’^ point in 
the Markov chain, a new point is determined by throwing 
a random number according to a probability density called 
proposal density. We use Gaussian proposal densities, where 
for each parameter the mean is given by the current param¬ 
eter value and the width is regularly adjusted as discussed 
below. 

The x^ for the (n -f 1 )'* point is then calculated and com¬ 
pared to the x^ for the point. If the new point shows 
better or equal agreement between predictions and measure¬ 
ments, 

xl+\ < Xl, ( 2 ) 

it is accepted. If the [n + 1)'^’ point shows worse agreement 
between the predictions and measurements, it is accepted 
with probability 


p = exp - 


yl _ y2 

Afi+l Ah 


(3) 


and rejected with probability 1 — p. If the (n + 1)'* point 
is rejected, new parameter values are generated based on the 
n'^ point again. If the (n + 1)^'' point is accepted, new param¬ 
eter values are generated based on the {n -f 1)'* point. Since 
the primary goal of using the MCMC method is the accurate 
determination of the best fit point and a high sampling den¬ 
sity around that point in the region of A x^ < 6, while allow¬ 
ing the MCMC method to escape from local minima in the 
X^ landscape, it is mandatory to neglect rejected points in 
the progression of the Markov chain. However, the rejected 
points may well be used in the frequentist interpretation of 
the Markov chain and for the determination of the p-value. 
Thus, we store them as well in order to increase the overall 
sampling density. 

An automated optimisation procedure was employed to 
determine the width of the Gaussian proposal densities for 


'Since the allowed region for all observable sets tested in Section]^ 
differ only marginally, it does not matter significantly which observable 
set is chosen for the initial scan, as long as it efficiently samples the 
relevant parameter space 
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each parameter for different targets of the acceptance rate of 
proposed points. Since the frequentist interpretation of the 
Markov chain does not make direct use of the point density, 
we can employ chains, where the proposal densities vary 
during their evolution and in different regions of the param¬ 
eter space. We update the widths of the proposal densities 
based on the variance of the last ^(500) accepted points in 
the Markov chain. Also, different ratios of proposal densities 
to the variance of accepted points are used for chains started 
in different parts of the parameter space, to optimally scan 
the widely different topologies of the surface at different 
SUSY mass scales. These differences stem from the varying 
degree of correlations between different parameters required 
to stay in agreement with the data, and from non-linearities 
between the parameters and observables. They are also the 
main reason for the excessive amount of points needed for a 
typical SUSY scan, as compared to more nicely behaved pa¬ 
rameter spaces. It has been ensured that a sufficient number 
of statistically independent chains yield similar scan results 
over the full parameter space. For the final interpretation, all 
statistically independent chains are added together. 

A total of 850 million valid points have been tested. The 
point with the lowest overall is identified as the 

best fit point. 

2.1.2 Interpretation of Markov chain results 

In addition to the determination of the best fit point it is also 
of interest to set limits in the cMSSM parameter space. For 
the Frequentist interpretation the measure 

^x^ = x^-xlm ( 4 ) 

is used to determine the regions of the parameter space 
which are excluded at various confidence levels. For this 
study the one dimensional 1(7 region (Ax^ < 1) and the two 
dimensional 2(7 region (Ax^ < 6) are used. In a Gaussian 
model, where all observables depend linearly on all parame¬ 
ters and where all uncertainties are Gaussian, this would cor¬ 
respond to the 1-dimensional 68% and 2-dimensional 95% 
Confidence Level (CL) regions. The level of observed de¬ 
viation from this pure Gaussian approximation shall be dis¬ 
cussed together with the results of the toy fits, which are an 
ideal tool to resolve these differences. 


2.2 Determining the p-value 


In all previous instances of SUSY fits, no true frequentist p- 
value for the fit is calculated. Instead, usually the x^,„/ndf is 
calculated, from which for a linear model with Gaussian ob¬ 
servables a p-value can easily derived. It has been observed 
that the Z^,>,/ndf of constrained SUSY model fits such as 
the cMSSM have been degrading while the direct limits 
on the sparticle mass scales from the LHC got stronger 


(see e.g. |3^^). Thus, there is the widespread opinion 


that the cMSSM is obsolete. However, as the cMSSM is 
a highly non-linear model and the observable set includes 
non-Gaussian observables, such as one-sided limits and the 
ATLAS 0-lepton search, it is not obvious that the Gaussian 
^^-distribution for ndf degrees of freedom can be used to 
calculate an accurate p-value for this model. Hence the main 
question in this paper is: How obsolete is the cMSSM, ex¬ 
actly? To answer this, a machinery to re-interpret the scan 
described above had to be developed, since re-scanning the 
parameter space for each individual toy observable set is 
computationally prohibitive at present. Because during this 
re-interpretation of the original scan a multitude of different 
cMSSM points might be chosen as optima of the toy fits, 
such a procedure sets high demands on the scan density also 
over the entire approximate 2 to 3 sigma region around the 
observed optimum. 

2.2.1 General Procedure 

After determining the parameter values that provide the best 
combined description of the observables suitable to con¬ 
strain the model, the question of the p-value for that model 
remains: Under the assumption that the tested model at the 
best fit point is the true description of nature, what is the 
probability p to get a minimum x^ as bad as, or worse than, 
the actual minimum 

For a set of observables with Gaussian uncertainties, 
this probability is calculated by means of the distri¬ 
bution and is given by the regularised Gamma function, 

P — P ^ ■ Here, n is the number of degrees of free¬ 

dom of the fit, which equals the number of observables mi¬ 
nus the number of free parameters of the model. 

In some cases, however, this function does not describe 
the true distribution of the x^- Reasons for a deviation in¬ 
clude non-linear relations between parameters and observ¬ 
ables (as evident in the cMSSM, where a strong variation of 
the observables with the parameters at low parameter scales 
is observed, while complete decoupling of the observables 
from the parameters occurs at high scales), non-Gaussian 
uncertainties as well as one-sided constraints, that in ad¬ 
dition might constrain the model only very weakly. Also, 
counting the number of relevant observables n might be non¬ 
trivial: For instance, after the discovery of the Higgs bo¬ 
son at the LHC, the limits on different Higgs masses set by 
the LEP experiments are expected to contribute only very 
weakly (if at all) to the total x^ in n fit of the cMSSM. This 
is because the measurements at the LHC indicate that the 
lightest Higgs Boson has a mass significantly higher than 
the lower mass limit set by LEP. In such a situation, it is not 
clear how much such a one-sided limit actually is expected 
to contribute to the distribution of x^ values. 

Eor the above reasons, the accurate determination of the 
p-values for the fits presented in this paper requires the con- 
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sideration of pseudo experiments or “toy observable sets”. 
Under the assumption that a particular best fit point provides 
an accurate description of nature, pseudo measurements are 
generated for each observable. Each pseudo measurement 
is based on the best fit prediction for the respective observ¬ 
able, taking into account both the absolute uncertainty on 
the best fit point, as well as the shape of the underlying 
probability density function. For one unique set of pseudo 
measurements, the fit is repeated, and a new best fit point is 
determined with a new minimum Xbf /■ 

This procedure is repeated ntoy times, and the number 
Hp of fits using pseudo measurements with ■ > Xbf 
determined. The p-value is then given by the fraction 

P=—■ (5) 

Wtoy 

This procedure requires a considerable amount of CPU time; 
the number of sets of pseudo measurements is thus limited 
and the resulting p-value is subject to a statistical uncer¬ 
tainty. Given the true p-value, 

Poo = Urn p, (6) 

ntoy->~ 

Hp varies according to a binomial distribution 
Z?(np|poo,ntoy), which in a rough approximation gives 
an uncertainty of 


Ap = 


p-il-p) 


^toy 

on the /7-value. 


(7) 


2.2.2 Generation of Pseudo Measurements for the cMSSM 


In the present fit of the cMSSM a few different classes of ob¬ 
servables have been used and the pseudo experiments have 
been generated accordingly. In this work we distinguish dif¬ 
ferent smearing procedures for the observables; 


a) For a Gaussian observable with best fit prediction 
and an absolute uncertainty at the best fit point, 
pseudo measurements have been generated by throwing 
a random number according to the probability density 
function 


P 



O 


V2 


— „p ■ exp 


•‘r-off 


( 8 ) 


b) For the measurements of the Higgs signal strengths and 
the Higgs mass, the smearing has been performed by 
means of the covariance matrix at the best fit point. 

c) For the ATLAS 0-Lepton search | |8T| (see Section [TT] ), 
the number of observed events has been smeared accord¬ 
ing to a Poisson distribution. The expectation value of 
the Poisson distribution has been generated for each toy 
by taking into account the nominal value and the sys¬ 
tematic uncertainty on both the background and signal 
expectation at the best fit point. The systematic uncer¬ 
tainties are assumed to be Gaussian. 


d) The best fit point for each set of observables features a 
lightest Higgs boson with a mass well above the LEP 
limit. Assuming the best fit point, the number of ex¬ 
pected Higgs events in the LEP searches is therefore 
negligible and has been ignored. For this reason, the 
LEP limit has been smeared directly assuming a Gaus¬ 
sian probability density function. 

2.2.3 Rerunning the Fit 

Due to the enormous amount of CPU time needed to accu¬ 
rately sample the parameter space of the cMSSM and calcu¬ 
late a set of predictions at each point, a complete resampling 
for each set of pseudo measurements is prohibitive. 

For this reason the pseudo fits have been performed us¬ 
ing only the points included in the original Markov chain, 
for which all necessary predictions have been calculated in 
the original scan. 

In addition, an upper cut on the (calculated with re¬ 
spect to the real measurements) of < 15 has been ap¬ 
plied to further reduce CPU time consumption. The cut is 
motivated by the fact, that in order to find a toy best fit 
point that far from the original minimum, the outcome of the 
pseudo measurements would have to be extremely unlikely. 
While this may potentially prevent a pseudo fit from finding 
the true minimum, tests with completely Gaussian toy mod¬ 
els have shown that the resulting distributions perfectly 
match the expected distribution for all tested numbers of 
degrees of freedom. 

As will be shown in Section |4.3| in general we observe 
a trend towards less pseudo data fits with high x^ values 
in the upper tail of the distribution than expected from the 
naive gaussian case. This further justifies that the Ax'^ < 15 
cannot be expected to bias the full result of the pseudo data 
fits. 

Nevertheless, the p-value calculated using the described 
procedure may be regarded as conservative in the sense that 
the true p-value may very well be even lower. Hence, if it is 
found below a certain threshold of e.g. 5%, it is not expected 
that there is a bias that the true p-value for infinite statistics 
is found at larger values. If for a particular toy fit the true 
best fit point is not included in the original Markov chain, 
the minimum x^ for that pseudo fit will be larger than the 
true minimum for that pseudo fit, which artificially increases 
the p-value. 


3 Observables 

The parameters of the cMSSM are constrained by preci¬ 
sion observables, like (g — 2)p, astrophysical observations 
including in particular direct dark matter detection limits 
and the dark matter relic density, by collider searches for 
supersymmetric particles and by the properties of the Higgs 
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Table 1 Precision observables used in the fit. 


- a™ 

(28.7±8.0) X 10-'“ 

(82l|8^ 

sin^ 9eff 

0.23113±0.00021 


84 


mt 

(173.34±0.27±0.71)GeV 


S5 


m^^ 

(80.385 ±0.015) GeV 


86 


Anis 

(17.719±0.036±0.023)ps-' 


S7 



(2.90 ± 0.70) X 10-5 




ilS(b i/) 

(3.43 ±0.21 ±0.07) X 10-^ 




-> tv) 

(1.05±0.25)x 10-“* 


S7 



Table 2 Standard Model parameters that have been fixed. Please note 
that irtf, and are MS masses at their respective mass scale, while for 
all other particles on-shell masses are used. 


l/Oem 128.952 

Gf (1.1663787 x IQ-^) GeV^^ 

tts 0.1184 

mz 91.1876 GeV 

m* 4.18 GeV 

1.77682 GeV 
ntc 1.275 GeV 


83 

12 

12 

12 

12 

12 

S7 


boson. In this section we describe the observables that enter 
our fits. The measurements are given in Section HD while 
the codes used to obtain the corresponding model predic¬ 
tions are described in Section|T2] 


3.1 Measurements and exclusion limits 


We employ the same set of precision observables as in our 
previous analysis Ref. but with updated measurements 
as listed in Tab.[2 They include the anomalous magnetic 
moment of the muon (g — 2)^ = a^, the effective weak mix¬ 
ing angle sin^ 0eff, the masses of the top quark and W boson, 
the Bs oscillation frequency A nig, as well as the branching 
ratios —>• /f/f), ^{B —)■ tv), and £§{b —)■ sy). The Stan¬ 

dard Model parameters that have been fixed are collected in 
Tab.|^ Note that the top quark mass m, is used both as an 
observable, as well as a floating parameter in the fit, since it 
has a significant correlation especially with the light Higgs 
boson mass. 

Dark matter is provided by the lightest supersymmetric 
particle, which we require to be solely made up of the neu- 
tralino. We use the dark matter relic density = 0.1187 ± 
0.0017 as obtained by the Planck collaboration | [90) and 
bounds on the spin-independent WIMP-nucleon scattering 


cross section as measured by the LUX experiment |911. 


Supersymmetric particles have been searched for at the 
LHC in a plethora of final states. In the cMSSM parame¬ 
ter region preferred by the precision observables listed in 
Tab. [T] the LHC jets plus missing transverse momentum 
searches provide the strongest constraints. We thus imple¬ 


ment the ATLAS analysis of Ref. |81| in our fit, as de 


scribed in some detail in 13^ . Furthermore we enforce 
the LEP bound on the chargino mass, m„± > 103.5 GeV 
192) . The constraints from Higgs searches at LEP are in¬ 
corporated through the extension provided by HlGGS- 
Bounds 4.1.1 | |93| - |9^ , which also provides limits on addi¬ 
tional heavier Higgs bosons. The signal rate and mass mea¬ 
surements of the experimentally established Higgs boson at 
125 GeV are included using the program HiGGSSIGNALS 
1.2.0 | [97l (see also Ref. | [98) and references therein). HlG- 
GSSlGNALS is a general tool which allows the test of any 
model with Higgs-like particles against the measurements 
at the LHC and the Tevatron. Therefore, its default observ¬ 
able set of Higgs rate measurements is very extensive. As 
is discussed in detail in Section [43] this provides maximal 
flexibility and sensitivity on the constraints of the allowed 
parameter ranges, but is not necessarily ideally tailored for 
goodness-of-fit tests. There, it is important to combine ob¬ 
servables which the model on theoretical grounds cannot 
vary independently. In order to take this effect into account, 
in our analysis we compare five different Higgs observable 
sets; 


Set 1 (Large Observable Set): This set is the default ob¬ 
servable set provided with HiGGSSiGNALS 1.2.0, con¬ 
taining in total 80 signal rate measurements obtained 
from the LHC and Tevatron experiments. It contains all 
available sub-channel / category measurements in the 
various Higgs decay modes investigated by the experi¬ 
ments. Hence, while this set is most appropriate for re¬ 
solving potential deviations in the Higgs boson coupling 
structure, it comes with a high level of redundancy. De¬ 
tailed information on the signal rate observables in this 
set can be found in Ref. 198) . Eurthermore, the set con¬ 
tains 4 mass measurements, performed by ATLAS and 
CMS in the h — 77 and h —>■ ZZ^*) 4 f channels. It 
is used as a cross-check for the derived observable sets 
described below. 

Set 2 (Medium Observable Set): This set contains 10 in¬ 
clusive rate measurements, performed in the channels 
h ^ YY, h ^ ZZ, h ^ WW, Vh V{bb) (V = W,Z), 
and h ^ Txhy ATLAS and CMS, listed in Tab. [^ As in 
Set 7, 4 Higgs mass measurements are included. 

Sets (Small Observable Set): In this set, the h —^ 77, 
h —^ ZZ and h WW channels are combined to 
a measurement of a universal signal rate, denoted 
h — > 77,ZZ,VTVT in the following. Together with the 
Vh ^ V (bb) and /z — > TT from Set 2, we have in total 
6 rate measurements. Eurthermore, in each LHC exper¬ 
iment the Higgs mass measurements are combined. The 
observables are listed in Tab. |4| 
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Table 3 Higgs boson mass and rate observables of Set 2 (Medium Ob- Table 4 Higgs boson mass and rate observables of Set 3 (Small Ob¬ 
servable Set). servable Set). 


Experiment, Channel 


observed p. observed m^, 


ATLAS,/!-> WW-> fvfv 100 
ATLAS, /i ZZ -> 4^ (Ito “ 
ATLAS, /i 77 ^ 

ATLAS, h -> TT (IM 


0.99tl]i> 

1 - 431 o :35 (124.3±Ll)GeV 

3-J^-0.28 
1 44 + 0.51 
3 -^^- 0.43 


(126.8 ± 0.9) GeV 


ATLAS, Vh -3 V{bb) 

104 

1 0.17+g 

CMS, h WW ->• Iviv 
CMS,/!-> ZZ-> 4f 101 

105 

0-7271];“ 

0.93+^ (125.6±0.6)GeV 


CMS, h 
CMS, h 
CMS, l//i->y(^||l06| 


■ 77 1102 
■ttITM 


o.78i|;:i? 

i-ooto:io 


(125.4±Ll)GeV 


Set 4 (Combined Observable Set): In this set we further 
reduce the number of Higgs observables by combin¬ 
ing the ATLAS and CMS measurements for the Higgs 
decays to electroweak vector bosons (V = VT,Z), pho¬ 
tons, Zt-quarks and T-leptons. These combinations are 
performed by fitting a universal rate scale factor jj. to 
the relevant observables from within Set 1. Furthermore, 
we perform a combined fit to the Higgs mass observ¬ 
ables of Set 1, yielding m/, = (125.73 ±0.45) GeV0The 
observables of this set are listed in Tab. |5] 

Set 5 (Higgs mass only): Here, we do not use any Higgs 
signal rate measurements. We only use one combined 
Higgs mass observable, which in our case is m/, = 
(125.73 ±0.45) GeV (see above). 


3.2 Model predictions 


We use the following public codes to calculate the predic¬ 
tions for the relevant observables: The spectrum is calcu¬ 
lated with SPheno 3.2.4 |108 109| . First the two-loop 
RGBs 1 110 1 are used to obtain the parameters at the scale 
Q — -^mj^ mj^. At this scale the complete one-loop correc¬ 
tions to all masses of supersymmetric particles are calcu¬ 
lated to obtain the on-shell masses from the underlying DR 
parameters |[TTT]. A measure of the theory uncertainty due 
to the missing higher-order corrections is given by varying 


^Note that the computing time needed for creating the pseudo-data fits 
presented in Section [43] means that the fits were starting to be per- 
formed significantly before the combined measurement of the Higgs 
boson mass = 125.09± 0.21 GeV by the ATLAS and CMS col¬ 
laborations was published |99| . We therefore performed our own com¬ 
bination, based on earlier results as published in |100j{102^ . Given the 
applied theory uncertainty on the Higgs mass prediction of = 

3 GeV, a shift of 0.64 GeV in the Higgs boson mass has a very small 
effect of Ax^ ~ ^(0.64^/3^) = 0.046, which is negligible in terms of 
the overall conclusions in this paper. 


Experiment, Channel observed p observed m;, 


ATLAS,/!->WW,ZZ, 77 100 
ATLAS, /! -> TT fl^ 

ATLAS, Vh ■ 


■ V{bb) |W| 


1 ' 50 + 0.21 
^••^^-0.18 
1 44 + 0.51 
^•^^-0.43 

A 17+0.67 

'- 0.63 


(125.5±0.8)GeV 


CMS,/!->lVlV,ZZ, 77 + O-SO+jf (125.7±0.6) GeV 


CMS, 106 

CMS, WiV(l>&) 

106 

0.78+“ 

LOO+iU 

' The combination of the CMS h —> tVVk.ZZ, 77 channels has been per 

formed with HiGGSSiGNALS using results from Ref. |101|105 

107 

The combined mass measurement for CMS is taken from Ref. 

102 


Table 5 Higgs boson mass and rate observables of Set 4 (Combined 
Observable Set). 


Experiment, Channel 

observed p 

observed m* 

ATLAS+CMS, h -> WW,ZZ 
ATLAS+CMS, /! -> 77 

o.94+:17 

1 - 1 + 0:20 

(125.73 ±0.45) GeV 

ATLAS+CMS, It -> TT 

i-ii+:i 

- 

ATLAS+CMS, Vh,tth -> bb 

0.69+:“ 

- 


the scale Q between (7/2 and 2Q. We find that the uncer¬ 
tainty on the strongly interacting particles is about 1-2%, 
whereas for the electroweakly interacting particles it is of 
order a few per mille p09| . 

Properties of the Higgs bosons as well as a^, 
Ants, sin^Peff and mw are obtained with FeynHiggs 
2.10.1 1 112) , which - compared to FeynHiggs 2.9.5 and 
earlier versions - contains a significantly improved calcu¬ 
lation of the Higgs boson mass |113| for the case of a 
heavy SUSY spectrum. This improves the theoretical uncer¬ 
tainty on the Higgs mass calculation from about 3-4 GeV in 
cMSSM scenarios | |114f|l 16l to about 2GeV p4l . 

The B-physics branching ratios are calculated by Su- 
PERlSO 3.3 un). We have checked that the results for the 
observables, that are also implemented in SPHENO agree 
within the theoretical uncertainties, see also |118[ for a com¬ 
parison with other codes. 

For the astrophysical observables we use Ml- 
CROMEGAs 3.6.9 ]119| |120[ to calculate the dark 
matter relic density and DarkSUSY 5.0.5 m via the 
interface program ASTROFiT |122[ for the direct detection 
cross section. 

For the calculation of the expected number of signal 
events in the ATLAS jets plus missing transverse momen¬ 
tum search, we use the Monte Carlo event generator Her- 
WIGH-H- p23) and a carefully tuned and validated version 
of the fast parametric detector simulation Delphes |124|. 
For tanj3 = 10 and Aq = 0, a fine grid has been produced 
in Mo and Mi/ 2 - In addition, several smaller, coarse grids 
have been defined in Aq and tan j3 for fixed values of Mq and 




































































Table 6 Theoretical uncertainties of the precision observables used in 
the fit. 



Mo (GeV) 


Fig. 1 Comparison of the emulation of the ATLAS 0-Lepton search 
with the published ATLAS result. In red dots we show the Atlas me¬ 
dian expected limit; the red lines denote the corresponding 1CT uncer¬ 
tainty. The central black line is the result of the FITTING implemen¬ 
tation described in the text. The upper and lower black curves are the 
corresponding 1CT uncertainty. The yellow dots are the observed Atlas 
limit. 


Mi /2 along the expected ATLAS exclusion curve to correct 
the signal expectation for arbitrary values of Aq and tan j3. 
We assume a systematic uncertainty of 10% on the expected 
number of signal events. In Fig. [T] we compare the expected 
and observed limit as published by the ATLAS collaboration 
to the result of our emulation. The hgure shows that the pro¬ 
cedure works reasonably well and is able to reproduce with 
sufficient precision the expected exclusion curve, including 
the ±1(7 variations. 

We reweight the events depending on their production 
channel according to NLO cross sections obtained from 
ProSPINO 0125H127| . Renormalisation and factorisation 
scales have been chosen such that the NLOh-NLL resummed 
cross section normalisations [128 ■ 132] are reproduced for 
squark and gluino production. 

For all predictions we take theoretical uncertainties into 
account, most of which are parameter dependent and reeval¬ 
uated at every point in the MCMC. For the precision ob¬ 
servables, they are given in Tab. For the dark matter relic 
density we assume a theoretical uncertainty of 10%, for the 
neutralino-nucleon cross section entering the direct detec¬ 
tion limits we assign 50% uncertainty (see Ref. for a 
discussion of this uncertainty arinsing from pion-nuclueon 
form factors), for the Higgs boson mass prediction 2.4%, 
and for Higgs rates we use the uncertainties given by the 
LHC Higgs Cross Section Working Group |133| . 

One common challenge for computing codes specifi¬ 
cally developed for SUSY predictions is that they might not 
always exactly predict the most precise predictions of the 
SM value in the decoupling limit. The reason is that spe- 
cihc loop corrections or renormalisation conventions are not 


Ou — 

1% 

sin Self 

0.05% 

m, 

IGeV 

mw 

0 .01% 

Am, 

24% 

SS{Bs 11 fi) 

26% 

S§{b —>■ jy) 

14% 

SsIb^tv) 

20 % 


always numerically implemented in the same way, or that 
SM loop contributions might be missing from the SUSY 
calculation. In most cases these differences are within the 
theory uncertainty, or can be used to estimate those. One 
such case of interest for this fit occurs in the program Feyn- 
Higgs, which does not exactly reproduce the SM Higgs de¬ 
coupling limit |134[ as used by the LHC Higgs cross-section 
working group To compensate this, we rescale the 

Higgs production cross sections and partial widths of the 
SM-like Higgs boson. We determine the scaling factors by 
the following procedure ]134) : We hx tan)3 = 10. We set 
all mass parameters of the MSSM (including the param¬ 
eters jj. and niA of the Higgs sector) to a common value 
'«SUSY- We require all sfermion mixing parameters Ay to be 
equal. We vary them by varying the ratio Yf/msusv, where 
X, =At — jJ ./tan j3. The mass of the Higgs boson becomes 
maximal for values of this ratio of about ±2. We scan the 
ratio between these values. In this way we hnd for each 
'WSUSY two parameter points which give a Higgs boson mass 
of about 125.5 GeV. One of these has negative X,, the other 
positive Xf. We then determine the scaling factor by requir¬ 
ing that for OTsusY = 4 TeV and negative Xt the production 
cross sections and partial widths of the SM-like Higgs bo¬ 
son are the same as for a SM Higgs boson with the same 
mass of 125.5 GeV. We then determine the uncertainty on 
this scale factor by comparing the result with scale factors 
which we would have gotten by choosing msusY = 3 TeV, 
»*SUSY = 5 TeV or a positive sign of X,. This additional un¬ 
certainty is taking into account in the computation. By 
this procedure we derive scale factors between 0.95 and 1.23 
with uncertainties of less than 0.6%. 


4 Results 


In Section |4T[ we show results based on the simplistic and 
common prohle likelihood technique, which all frequentist 
fits, including us, have hitherto been employing. In Sec¬ 
tion |4.2| a full scan of the allowed parameter space for a 
stable vacuum is shown, before moving on to novel results 


from toy fits in Section 4.3 
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(a) 1 O' and 2a contour plot in in the (Mo,Mi/2)“Pl4ne for the 
Small Observable Set. 



0^^2000^^4CjW^^60W^^8000^^10^ 

HA, (GeV) 


(b) 1 a and 2 a contour plot in in the (Mo,Mi/ 2 )-plane for the 
Medium Observable Set. 



0 2000 4000 6000 8000 10000 

M„ (GeV) 


(c) 1 a and 2a contour plot in the (Mo,Mi/ 2 )-plane for the 
Large Observable Set. 



(d) 1 (7 and 2 a contour plot in the (Aq , tan /3 )-plane for the 
Medium Observable Set. 


Fig. 2 1 a and 2 a contour plots for different projections and different observable sets. It can be seen that the preferred parameter region does not 
depend on the specific observable set. 


4.1 Profile likelihood based results 


In this section we describe the preferred parameter space re¬ 
gion of the cMSSM and its physical properties. Since a truly 
complete frequentist determination of a conhdence region 
would require not only to perform toy hts around the best 
ht point (as described in Section 2.2 and |4.3| l but around ev¬ 
ery cMSSM parameter point in the scan, we rely here on 
the prohle likelihood technique. This means, we show vari¬ 
ous projections of the lD-l(j/lD-2ff/2D-2a regions dehned 
as regions which satisfy Ax^ < 1/4/5.99 respectively. In 
this context, prohle likelihood means that out of the 5 phys¬ 
ical parameters in the scan, the parameters not shown in a 
plot are prohled, i.e. a scan over the hidden dimensions is 
performed and for each selected visible parameter point the 
lowest x'^ value in the hidden dimensions is chosen. Ob¬ 
viously, no systematics nuisance parameters are involved, 
since all systematic uncertainties are given by relative or 
absolute Gaussian uncertainties, as discussed in Section]^ 
One should keep in mind that this correspondence is actu¬ 


ally only exact when the observed distribution of values 
in a set of toy hts is truly ^^-distributed, which - as dis¬ 
cussed in Section |43] - is not the case. Nevertheless, since 
the exact method is not computationally feasible, this stan¬ 
dard method, as used in the literature in all previous frequen¬ 
tist results, gives a reasonable estimation of the allowed pa¬ 
rameter space. In Section [43| more comparisons between the 
sets of toy ht results and the prohle likelihood result will be 
discussed. 

Note that for the discussion in this and the next section, 
we treat the region around the best ht point as “allowed”, 
even though, depending on the observable set, an exclusion 
of the complete model will be discussed in Section |43] 

All hve Higgs input parameterisations introduced in 
Section]^ lead to very similar results when interpreted with 
the prohle likelihood technique. As an example, Figures [2(aj| 
-|2(c)|show the (Mq, M]/ 2 )-projection of the best ht point, 
1D-1(7 and 2D-2(7 regions for the Small, the Large and the 
Medium Observable Set. Thus, in the remainder of this sec- 
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Table 7 Central values and 1 a uncertainties of the free model param¬ 
eters at the global and secondary minimum when using the Medium 
Observable Set 


Parameter 

global minimum 

secondary minimum 

Mo (GeV) 

387.4+t8|;7 

8983.4i?*,s, 

Mi /2 (GeV) 

918.2+1” 

2701 . 1 ”^“ 

Ao (GeV) 

-2002.8”^”^ 

5319.0+™9| 

tan/I 

17 . 7 ^ 1^1 


m, (GeV) 

174.3+j;[ 

172. 


> 




8000 10000 
Mo (GeV) 


Fig. 3 2 a region in the (Mq, Mj/ 2 )-plane for the Medium Observable 
Set. Regions with different dark matter annihilation mechanisms are 
indicated. The enclosed red areas denote the best fit regions shown in 
Fig.|2(b)l 


tion, we concentrate on results from the Medium Observable 
Set. 


The (Mo, Mi/ 2 )-projection in Fig. 2(b) 


shows two dis¬ 
joint regions. While in the region of the global minimum, 
values of less than 900 GeV for Mo and less than 1300 GeV 
for Mj 12 are preferred at 1 a, in the region of the secondary 
minimum values of more than 7900 GeV for Mo and more 
than 2100 GeV forMi /2 are favored (see also Tab.j^. 

The different regions are characterised by different dom¬ 
inant dark matter annihilation mechanisms as shown in 
Fig. Here we define the different regions similarly to 
Ref. (4§ by the following kinematical conditions, which we 
have adapted such that each point of the 2D-2(J region be¬ 
longs to at least one region; 

- fi coannihilation: mf, /m„o — 1 <0.15 

‘ /Cl 

- fi coannihilation: mj — 1 < 0.2 

- Xi coannihilation; m^± ~ 1 <0.1 


-A///funnel; |myi/2m„o — 11 < 0.2 

A1 

- focus point region; |/f /m^o — 11 < 0.4 

Z[ 

With these definitions each parameter point of the pre¬ 
ferred 2D-2(7 region belongs either to the fi coannihilation 


or the focus point region. Additionally a subset of the points 
in the fi coannihilation region fulfills the criterion for the 
A/M-funnel, while some points of the focus point region 
fulfill the criterion for xf" coannihilation. No point in the 
preferred 2D-2a region fulfills the criterion for ti coannihi¬ 
lation, due to relatively large stop masses. 

At large Mq and low Mj /2 a thin strip of our preferred 
2D-2(7 region is excluded at 95% confidence level by AT¬ 
LAS jets plus missing transverse momentum searches re¬ 
quiring exactly one lepton |135| or at most one lepton and 
b-jets ]136| in the final state. Therefore an inclusion of these 
results in the fit is expected to remove this small part of the 
focus point region without changing any conclusion. 

Also note that the parameter space for values of Mq 
larger than lOTeV was not scanned such that the preferred 
20-2(7 focus point region is cut off at this value. Because 
the decoupling limit has already been reached at these large 
mass scales we do not expect significant changes in the pre¬ 
dicted observable values when going to larger values of Mq. 
Hence we also expect the ID-Iff region to extend to larger 
values of Mq than visible in Fig. |2(b) due to a low sampling 
density directly at the 10 TeV boundary. For the same reason 
this cut is not expected to influence the result of the p-value 
calculation. If it does it would only lead to an overestimation 
of the /7-value. 

In the fi coannihilation region negative values of Aq be¬ 
tween —4000 GeV and —1400 GeV and moderate values of 
tanjS between 6 and 35 are preferred, while in the focus 
point region large positive values of Aq above 3400 GeV 
and large values of tan j3 above 36 are favored. This can be 
seen in the (Aq, tanj3)-projection shown in Fig. 2(d) and in 
Tab. El 

While the fi coannihilation region predicts a spin inde¬ 
pendent dark matter-nucleon scattering cross section which 
is well below the limit set by the LUX experiment, this mea¬ 
surement has a significant impact on parts of the focus point 
region for lightest supersymmetric particle (LSP) masses be¬ 
tween 200 GeV and 1 TeV, as shown in Fig.j^ The plot also 
shows how the additional uncertainty of 50% on ffsi shifts 
the implemented limit compared to the original limit set by 
LUX. It can be seen that future improvements by about 2 
orders of magnitude in the sensitivity of direct detection ex¬ 
periments, as envisaged e.g. for the future of the XENON IT 
experiment 1 137| , could at least significantly reduce the re¬ 
maining allowed parameter space even taking the systematic 
uncertainty into account, or finally discover SUSY dark mat¬ 
ter. 


The predicted mass spectrum of the Higgs bosons and 
supersymmetric particles at the best fit point and in the one¬ 
dimensional Iff and 2(7 regions is shown in Fig. Due to 
the relatively shallow minima of the fit a wide ranges of 
masses is allowed at 2ff for most of the particles. The masses 
of the coloured superpartners are predicted to lie above 1.5 
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Fig. 4 IDlff and 2D2(7 regions in m^o-asi for the Medium Observ- 

A.\ 

able Set. The LUX exclusion is shown in addition. The dashed line 
indicates a contribution of 1.64 which corresponds to a 90% CL 
upper limit. This line does not match the LUX exclusion line because 
we use a theory uncertainty of 50% on Osi- 



118 120 122 124 126 128 130 
m, (GeV) 



hO ^0 H" H" x“ X" Z" X“ x; X^ Tr I X, \ b, t, t^ 9 


Fig. 5 The Higgs and supersymmetric particle mass spectrum as pre¬ 
dicted by our fit using the Medium set of Higgs observables. 


Fig. 6 Our predicted mass of the light Higgs boson, together with the 
1 (7 and 2 (T ranges. The LHC measurements used in the fit are shown as 
well. Note that the correlated theory uncertainty of = 3 GeV is 

not shown in the the plot. The relative smallness of the 68% CL region 
of the fitted mass of = 1.1 GeV is caused by constraints from 

other observables. 


diets a slightly smaller cross section in all channels except 
the bbh channel. The predicted signal strengths fi in the dif¬ 
ferent final states for p-p collisions at a centre-of-mass en¬ 
ergy of 8 TeV is also slightly smaller than the SM prediction, 
as shown in Fig. The current precision of these measure¬ 
ments does, however, not allow for a discrimination between 
the SM and the cMSSM based solely on measurements of 
Higgs boson properties. 


4.2 Vacuum stability 


TeV, but due to the focus point region also masses above 10 
TeV are allowed at 1(7. Similarly the heavy Higgs bosons 
have masses of about 1.5 TeV at the best fit point, but masses 
of about 6 TeV are preferi'ed in the focus point region. The 
sleptons, neutralinos and charginos on the other hand can 
still have masses of a few hundred GeV. 

A lightest Higgs boson with a mass as measured by the 
ATLAS and CMS collaborations can easily be accommo¬ 
dated, as shown in Fig.|^ As required by the signal strength 
measurements, it is predicted to be SM-like. Fig.|7]shows a 
comparison of the Higgs production cross sections for dif¬ 
ferent production mechanisms in p-p collisions at a centre- 
of-mass energy of 14 TeV. These contain gluon-fusion (ggh), 
vector boson fusion (qqh), associated production (Wh, Zh), 
and production in assiciation with heavy quark flavours 
(bbh, tth). Compared to the SM prediction, the cMSSM pre- 


The scalar sector of the SM consists of just one complex 
Higgs doublet. In the cMSSM the scalar sector is dramat¬ 
ically expanded with an extra complex Higgs doublet, as 
well as the sfermions clk, Vl, ul,Rj d^R of the first family, 
and correspondingly of the second and third families. Thus 
there are 25 complex scalar fields. The corresponding com¬ 
plete scalar potential Vcmssm is fixed by the five parameters; 
(Mo, Mi/ 2 , A, tanjS, sgn(/i)). The minimal potential energy 
of the vacuum is obtained for constant scalar field values 
everywhere. Given a fixed set of these cMSSM parameters, 
it is a computational question to determine the minimum of 
ycMSSM- Ideally this minimum should lead to a Higgs vac¬ 
uum expectation such that SU(2)iXU(l)y —>'U(1 )em, as in 
the SM. However, it was observed early on in supersym¬ 
metric model building, that due to the extended scalar sec¬ 
tor, some sfermions could obtain non-vanishing vacuum ex- 
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I 68 % CL 
95 % CL 


— best fit value 
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'mssf 



a/a 


,SM , 


14TeV 


Fig. 7 Predicted production cross sections at 14 TeV of the light Higgs 
boson relative to the SM value for a Higgs boson with the same mass 


pectation values (vevs). There could be additional minima 
of the scalar potential which would break SU(3)c and/or 
Lf(l)EM and thus colour and/or charge |7l |138 - [140) . If these 
minima are energetically higher than the conventional elec- 
troweak breaking minimum, then the latter is considered 
stable. If any of these minima are lower than the conven¬ 
tional minimum, our universe could tunnel into them. If 
the tunneling time is longer than the age of the Universe 
of 13.8 gigayears |j^, we denote our favored vacuum as 
metastable, otherwise it is unstable. However, this is only 
a rough categorisation. Since even if the tunneling time is 
shorter than the age of the universe, there is a finite proba¬ 
bility, that it will have survived until today. When computing 
this probability, we set a limit of 10% survival probability. 
We wish to explore here the vacuum stability of the preferred 
parameter ranges of our fits. 

The recent observation of the large Higgs boson mass 
requires within the cMSSM large stop masses and/or a 
large stop mass splitting. Together with the tuning of the 
lighter stau mass to favor the stau co-annihilation region 
(for the low Mq fit region), this typically drives fits to fa¬ 
vor a very large value of |Ao| relative to \Mq\, cf Tab. |7] 
(For alternative non-cMSSM models with a modified stop 
sector, see for example 1141 ■ 144|.) This is exactly the re¬ 


1 68 % CL 
95 % CL 


— best fit value 
-^data 




ATL, h ^ WW ^ Iviv 
ATL, h ^ ZZ ^ 41 
ATL, h ^ YY 
ATL, h ^ TT 
ATL, Vh ^ Vbb 
CMS, h ^ WW ^ Iviv 
CMS, h ^ ZZ ^ 41 
CMS, h ^ yY 
CMS, h ^ XT 
CMS, Vh ^ Vbb 




Fig. 8 Our predicted fl values of the light Higgs boson relative to the 
SM value for a Higgs boson with the same mass. The measurements 
used in the fit are shown as well. 


gion, which typically suffers from the SM-like vacuum be¬ 
ing only metastable, decaying to a charge- and/or colour¬ 
breaking (CCB) minimum of the potential 1145 146) . 

For the purpose of a fit, in principle a likelihood value 
for the compatibility of the lifetime of the SM-like vacuum 
of a particular parameter point with the observation of the 
age of the Universe should be calculated and should be im¬ 
plemented as a one-sided limit. Unfortunately, the effort re¬ 
quired to compute all the minima of the full scalar potential 
and to compute the decay rates for every point in the MCMC 
and to implement this in the likelihood function is beyond 


present capabilities 11451. 


Effectively, whether or not a parameter point has an un¬ 
acceptably short lifetime has a binary yes/no answer. There¬ 
fore, as a first step, and in the light of the results of the possi¬ 
ble exclusion of the cMSSM in Section [43j we overlay our 
fit result from Section ItTI over a scan of the lifetime of the 
cMSSM vacuum over the complete parameter space. 

The systematic analysis of whether a potential has min¬ 
ima which are deeper than the desired vacuum configura¬ 


tion has been automated in the program Vevacious [ 1471. 
When restricting the analysis to only a most likely relevant 
subset of the scalar fields of the potential, i.e. not the full 
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Fig. 9 Preferred ID-lff and 2D-2c7 regions in M 0 -M 1/2 for the 
Medium Observable Set. The filled areas contain stable points, while 
the doted lines enclose points which are metastable but still might be 
very long-lived. The whole preferred 2D-2a focus point region leads 
to a stable vacuum, while the coannihilation region contains both sta¬ 
ble and metastable points. There are no stable points in the preferred 
ID-1(7 coannihilation region. 


25 complex scalar fields, and ignoring the calculation of 
lifetimes, this code runs sufficiently fast that we are able 
to present an overlay of which parameter points have CCB 
minima deeper than the SM-like minimum in Figures]^ and 
10 However, only the stop and stau fields were allowed to 


have non-zero values in determining the overlays, in addi¬ 
tion to the neutral components of the two complex scalar 
Higgs doublets. The Tl r, Ilr are suspected to have the 
largest effect 11451. The computation time when including 
more scalar fields which are allowed to vary increases expo¬ 
nentially. Thus the more detailed investigations below are re¬ 
stricted to a set of benchmark points. Note that the overlays 
in Figures]^ and 10 only show whether metastable vacua 
might occur at a given point, or whether the vacuum is in¬ 
stable at all. The actual lifetime is not yet considered in this 
step. See the further considerations below. 

There are analytical conditions in the literature for 
whether MSSM parameter points could have dangerously 
deep CCB minima, see for example |7 |139 140[ 148 ■ 1521. 
These can be checked numerically in a negligible amount of 
CPU time, while performing a fit. However, these conditions 
have been explicitly shown to be neither necessary nor suffi¬ 
cient 11531. In particular they have also been shown numer¬ 
ically to be neither necessary nor sufficient for the relevant 
regions of the cMSSM parameter space which we consider 
here |145| . 

Since the exact calculation of the lifetime of a metastable 
SM-like vacuum is so computationally intensive, we unfor¬ 
tunately must restrict this to just the best-fit points of the 
stau co-annihilation and focus point regions of our the fit, as 
determined in Section 4.1 As an indicator, though, the ex¬ 
tended fi co-annihilation region of the cMSSM investigated 


Fig. 10 Preferred ID-lcJ and 2D-2(T regions in tan/l-Ao/Mo for the 
Medium Observable Set. The filled areas contain stable points, while 
the doted lines enclose points which are metastable but still might be 
very long-lived. Points leading to a metastable vacuum have usually 
larger negative values of Ao relative to Mq when compared to points 
with a stable vacuum at the same tan/I. The part of the ID-1(7 region 
which belongs to the focus point region fulfills Aq/Mq 0 and is stable, 
while the part which belongs to the coannihilation region consists of 
points with relatively large negative values of Aq /Mq and is metastable. 


in GH) had SM-like vacuum lifetimes, which were all ac¬ 
ceptably long compared to the observed age of the Universe. 


The ID 1(7 best-fit points in Section 4.1 where checked 
for undesired minima, allowing, but not requiring, simulta¬ 
neously for all the following scalar fields to have non-zero, 
real VEVs: H^,H^,fL,fR,v^L,bL,bRjLjR- The focus point 


region best-fit point was found to have an absolutely stable 
SM-like minimum against tunneling to other minima, as no 
deeper minimum of Vcmssm was found at the 1 -loop level. 
The SM-like vacuum of the best-fit fi co-annihilation point 
was found to be metastable, with a deep CCB minimum 
with non-zero stau and stop VEVs. Eurthermore there were 
unbounded-from-below directions with non-zero values for 
the /r-sneutrino scalar field in combination with nonzero 
values for both staus, or both sbottoms, or both chiralities 
of both staus and sbottoms. This does not bode very well 
for the absolute best-fit point of our cMSSM fit. However, 
further effects must be considered. 


The parameter space of the MSSM which has directions 
in field space, where the tree-level potential is unbounded 
from below was systematically investigated in Ref. |152|. 
We confirmed the persistence of the runaway directions 
at one loop with Vevacious out to field values of the 
order of twenty times the renormalisation scale. This is 
about the limit of trustworthiness of a fixed-order, fixed- 
renormalisation-scale calculation 047). However, this is not 
quite as alarming as it may seem. The appropriate renor¬ 
malisation scale for very large field values should be of the 
order of the field values themselves, and for field values of 
the order of the GUT scale, the cMSSM soft SUSY-breaking 
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mass-squared parameters by definition are positive. Thus the 
potential at the GUT scale is bounded from below, as none of 
the conditions for unbounded-from-below directions given 
in m can be satisfied without at least one negative mass- 
squared parameter. Note, even the Standard Model suffers 
from a potential which is unbounded from below at a fixed 
renormalisation scale. Though in the case of the SM it only 
appears at the one-loop level. Nevertheless, RGBs show that 
the SM is bounded from below at high energies |T54). 

Furthermore, the calculation of a tunneling time out of 
a false minimum does not technically require that the Uni¬ 
verse tunnels into a deeper minimum. In fact, the state which 
dominates tunneling is always a vacuum bubble, with a field 
configuration inside, which classically evolves to the true 
vacuum after quantum tunneling 1 155[ 156) . Hence the life¬ 
time of the SM-like vacuum of the stau co-annihilation best- 
fit point could be calculated at one loop even though the 
potential is unbounded from below at this level. The min¬ 
imal energy barrier through which the SM-like vacuum of 
this point can tunnel is associated with a final state with 
non-zero values for the scalar fields i/®, f/., Xr, and 

The lifetime was calculated by using the program Veva- 


CIOUS through the program CosmoTransitions 11571 to 
be roughly g^ooo ^ times the age of the Universe. 

Therefore, we consider the X\ coannihilation region best-fit 
point as effectively stable. 

As well as asking whether a metastable vacuum has a 
lifetime at least as long as the age of the Universe at zero 
temperature, one can also ask whether the false vacuum 
would survive a high-temperature period in the early Uni¬ 
verse. Such a calculation has been incorporated into Veva- 


CIOUS 11581. In addition to the fact that the running of the 
Lagrangian parameters ensures that the potential is bounded 
from below at the GUT scale, the effects of non-zero tem¬ 
perature serve to bound the potential from below, as well. In 
fact the CCB minima of Vcmssm evaluated at the parameters 
of the stau co-annihilation best-fit point are no longer deeper 
than the configuration with all zero VEVs, which is assumed 
to evolve into the SM-like minimum as the Universe cools, 
for temperatures over about 2300 GeV. The probability of 
tunneling into the CCB state integrated over temperatures 
from 2300 GeV down to 0 GeV was calculated to be roughly 
exp(e^^*'‘'*’). So while having a non-zero-temperature de¬ 
cay width about ^-2000^^-4000 _ ^+2000 ^jj^es larger than 
the zero-temperature decay width, the SM-like vacuum, or 
its high-temperature precursor, of the stau co-annihilation 
best-fit point has a decay probability which is still utterly 
insignificant. 


4.3 Toy based results 

Pseudo datasets have been generated for a total of 7 differ¬ 
ent minima based on 6 different observable sets. For the 


Medium, Small and Combined Observable Sets, roughly 
1000 sets of pseudo measurements have been taken into ac¬ 
count, as well as for the observable set without the Higgs 
rates. For the Medium Observable Set, in addition to the 
best fit point, we also study the p-value of the local mini¬ 
mum in the focus point region. Due to relaxed requirements 
on the statistical uncertainty of a p-value in the range of 
^(0.5), as compared to ^(0.05), we use only 125 pseudo 
datasets for the Large Observable Set. Finally, to study the 
importance of a total of 500 pseudo datasets have 

been generated based on the best fit point for the Medium 
Observable Set without A summary of all p-values 

with their statistical uncertainties and a comparison to the 
naive p-value according to the ^^-distribution for Gaussian 
distributed variables is shown in Tab. 0 

Figure 11(a) shows the ^^-distribution for the best fit 


point of the Medium Observable Set, from which we de¬ 
rive a p-value of (4.9 ± 0.7)%. As a comparison we also 
show the ^^-distribution for the pseudo fits using the Com¬ 
bined Observable Set in figure |1 l(b)| Both distributions 
are significantly shifted towards smaller j^-values com¬ 
pared to the corresponding ^^-distributions for Gaussian 
distributed variables. Several observables are responsible for 
the large deviation between the two distributions, as shown 


in Fig. 12(a) where the individual contributions of all ob¬ 
servables to the minimum of all pseudo best fit points are 
plotted. 

First, HiggsBounds does not contribute significantly 
to the at any of the pseudo best-fit points, which is also 
the case for the original fit. The reason for this is, that for the 
majority of tested points, the x^ contribution from HlGGS- 
Bounds reflects the amount of violation of the LEP limit 
on the lightest Higgs boson mass by the model. Since the 
measurements of the Higgs mass by ATLAS and CMS lie 
significantly above this limit, it is extremely unlikely that 
in one of the pseudo datasets the Higgs mass is rolled such 
that the best fit point would receive a penalty due to the 
LEP limit. This effectively eliminates one degree of free¬ 
dom from the fit. In addition, the predicted masses of A,H 
and lie in the decoupled regime of the allowed cMSSM 
parameter space. Thus there are no contributions from heavy 
Higgs or charged Higgs searches as implemented in HlGGS- 
Bounds. 

The same effect is observed slightly less pronounced for 
the LHC and LUX limits, where the best fit points are much 
closer to the respective limits than in the case of HlGGS- 
Bounds. Finally we observe that for each pseudo dataset 
the cMSSM can very well describe the pseudo measurement 
of the dark matter relic density, which further reduced the 
effective number of degrees of freedom. 


Figure 12(a) also shows that the level of disagreement 
between measurement and prediction for (g — 2)^ is smaller 
in each single pseudo dataset than in the original fit with the 
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Table 8 Summary of p-values 


Observable Set 

X^/ndf 

naive p-value (%) 

toy p-value (%) 

Small 

27.1/16 

4.0 

1.9±0.4 

Medium 

30.4/22 

10.8 

4.9±0.7 

Combined 

17.5/13 

17.7 

8.3±0.8 

Medium (Focus Point) 

30.8/22 

10.0 

7.8±0.8 

Medium without (g-2) 

18.1/21 

64.1 

51±3 

Observable Set without Higgs rates 

15.5/9 

7.8 

1.3±0.4 



(a) Minimal values from toy fits using the 
Medium Observable Set. 



(b) Minimal values from toy fits using the 
Combined Observable Set. 


Fig. 11 Distribution of minimal values from toy fits using two different sets of Higgs observables. A distribution for Gaussian distributed 
variables is shown for comparison. 


real dataset. The 1-dimensional distribution of the pseudo 
best fit values of (g — 2)^ is shown in Fig. 13(a) The figure 
shows that under the assumption of our best fit point, not a 
single pseudo dataset would yield a prediction of {g — 2)^ 
that is consistent with the actual measurement. As a com¬ 


parison, Fig. 13(b) shows the 1-dimensional distribution for 
the dark matter relic density, where the actual measurement 
can well be accommodated in any of the pseudo best fit sce¬ 
narios. To further study the impact of (g — 2)^ on the p- 
value, we repeat the toy fits without this observable and get 
a p-value of (51 ± 3)%. This shows that the relatively low 
p-value for our baseline fit is mainly due to the incompatibil¬ 
ity of the (g — 2)^ measurement with large sparticle masses, 
which are however required by the LHC results. 

Interestingly, under the assumption that the minimum in 
the focus point region is the true description of nature, we 
get a slightly better p-value (7.8%) than we get with the ac¬ 
tual best fit point. Figure 12(b)| shows the individual contri¬ 
butions to the pseudo best fit at the pseudo best fit points 
for the toy fits performed around the local minimum in the 
focus point region. There are two variables with higher av¬ 
erage contributions compared to the global minimum: irifop 
and the LHC SUSY search. In particular for the LHC SUSY 
search, the LHC contribution to the total is, on average, 
significantly higher than for the pseudo best fit points for the 
global minimum. The number of expected signal events for 


the minimum in the focus point region is 0, while it is > 0 
for the global minimum. Pseudo best fit points with smaller 
values of the mass parameters, in particular pseudo best fit 
points in the f-coannihilation region, tend to predict an ex¬ 
pected number of signal events larger than zero. Since for 
the pseudo measurements based on the minimum in the fo¬ 
cus point region an expectation of 0 is assumed, this natu¬ 
rally leads to a larger contribution from the ATLAS 0- 
£ analysis. The effect on the distribution of the total is 
shown in figure Another reason might be that the fo¬ 
cus point region is sampled more coarsely than the region 
around the global minimum. This increases the probability 
that the fit of the pseudo dataset misses the actual best fit 
point, due to our procedure of using only the points in the 
original MCMC. This effect should however only play a mi¬ 
nor role, since the parameter space is still finely scanned and 
only a negligible fraction of scan points are chosen numer¬ 
ous times as best fit points in the pseudo data fits. 

To further verify that this effect is not only caused by the 
coarser sampling in the focus point region, we performed an¬ 
other set of 500 pseudo fits based on the global minimum, re¬ 
ducing the point density in the f-coannihilation region such 
that it corresponds to the point density around the local min¬ 
imum in the focus point region. We find that the resulting 
distribution is slightly shifted with respect to the x^ distri¬ 
bution we get from the full MCMC. The shift is, however. 














































16 
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(a) Toy fits smeared around the global minimum in the 
•f-coannihilation region. 


(b) Toy fits smeared around the local minimum in the focus point 
region. 


Fig. 12 Individual contributions of all observables/observable sets at the best fit points of the toy fits using the Medium Set of Higgs observables 
with observables smeared around the global and the local minimum of the observed contour. The predicted measurements at the best fit points of 
the individual pseudo data fits are used to derive the local CL intervals shown in the plots. These are compared with the individual contribution 
of each observable at the global or local minimum. Note the different scale shown on the top which is used for HiggsSignals, which contains 14 
observables. Also note that m* contains contributions from 4 measurements for this observable set. 
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(a) Anomalous magnetic moment of the muon. 


(b) Dark matter relic density. 


Fig. 13 Distribution of the predictions of the best fit points of the pseudo data fits for two different observables used in the fit, compai'ed with the 
respective measurements. 


too small to explain the difference between p-values we find 
for the global minimum and the local minimum in the focus 
point region. 

As an additional test, we investigate a simple toy model 
with only Gaussian observables and a single one-sided limit 
corresponding to the LHC SUSY search we use in our fit 
of the cMSSM. Also in this very simple model we find sig¬ 
nificantly different distributions for fits based on points 
in a region with/without a significant signal expectation for 
the counting experiment. We thus conclude that the true p- 
value for the local minimum in the focus point region is in 
fact higher than the true p-value for the global minimum of 
our fit. 

In order to ensure that there are no more points with a 
higher and a higher p-value than the local minimum in 


the focus point region, we generate 200 pseudo datasets for 
two more points in the focus point region. The two points 
are the points with the highest/lowest Mq in the local la 
region around the focus point minimum. We find that the 
distributions we get from these pseudo datasets are in good 
agreement with each other and also with the y} distribution 
derived from the pseudo experiments around the focus point 
minimum, and hence conclude that the local minimum in 
the focus point region is the point with the highest p-value 
in the cMSSM. 

To study the impact of the Higgs rates on the p-value, 
and in order to compare to the observable sets used by other 
fit collaborations, which exclude the Higgs rate measure¬ 
ments from the fit on the basis that in the decoupling regime 
they do not play a vital role, we perform toy fits for the 
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Fig. 14 Comparison of the x,'^ distributions obtained from toy fits using 
the global minimum and the local minimum in the focus point region 
of the Medium Observable Set. 


observable set without Higgs rates and derive a p-value of 
(1.3 ±0.4)%. In the decoupling limit, the cMSSM predic¬ 
tions for the Higgs rates are very close to the SM, so that 
the LHC is not able to distinguish between the two models 
based on Higgs rates measurements (see Fig. [^. Because 
of the overall good agreement between the Higgs rate mea¬ 
surements and the SM prediction, the inclusion of the Higgs 
rates in the fit improves the fit quality despite some tension 
between the ATLAS and CMS measurements. 


As discussed in Section]^ it is important to understand 
the impact of the parametrisation of the measurements on 
the p-value. To do so, we compare our baseline fit with two 
more extreme choices. First, we use the Small Observable 
Set which combines h ^ yy, h ^ ZZ, and h —> WW mea¬ 
surements but keeps ATLAS and CMS measurements sepa¬ 
rately. We use this choice because an official ATLAS combi¬ 
nation is available. The equivalent corresponding CMS com¬ 
bination is produced independently by us. Using this observ¬ 
able set we get a p-value of (1.9 ± 0.4)%. Here the cMSSM 
receives a penalty from the trend of the ATLAS signal 
strength measurements to values /i > 1 and of the CMS mea¬ 
surements towards /i < 1 in the three h —^ VV channels. 


As a cross-check, we employ the Large Observable Set, 
which contains all available sub-channel measurements sep¬ 
arately. Using this observable set, we get a p-value from 
the pseudo data fits of (41.6 ±4.4)%. As observed in Sec¬ 
tion 4.1 the Large Observable Set yields the same preferred 


parameter region as the Small, Medium and Combined Ob¬ 
servable sets. Yet, its p-value from the pseudo data hts sig¬ 
nificantly differs. 


To explain this interesting result we consider a simplified 
example: For / = 1,..., A, let Xj be Gaussian measurements 
with uncertainties (7, and corresponding model predictions 
a,(P) for a given parameter point P. We assume that the 
measurements from x„ to xn are uncombined measurements 


of the same observable; then a, = for all i > n. There are 
now two obvious possibilities to compare measurements and 
predictions: 

- We can compare each of the individual measurements 
with the corresponding model predictions by calculating 

N /.. „ \ 2 

Zspiit = E 

i=\ 

This would correspond to an approach where the model 
is confronted with all avaialable observables, irrespec¬ 
tive of the question whether they measure independent 
quantities in the model or not. One example for such a 
situation would be the Large Observable Set of Higgs 
signal strength measurements, where several observ¬ 
ables measure different detector effects, but the same 
physics. 

- We can first combine the measurements xt, i > n to a 
measurement x which minimises the function 


N 

f{x) = E 


Xj-X 

Oi 


(9) 


and has an uncertainty of ff and then use this combina¬ 
tion to calculate 

2 

( 10 ) 


Zcombined 

1=1 


n— 1 / \ 2 

^ Xi — Qi 




This situation now corresponds to hrst calculating one 
physically meaningful quantity (e.g. a common signal 
strength for h^yy'm all VBF categories, and all gg^h 
categories) and only then to confront the model to the 
combined measurement. 

Plugging in the explicit expressions for (x, a), using l/a^ = 
and defining = /W finds 

^combined Xsplit Xdata- (11) 

Hence doing the combination of the measurements be¬ 
fore the ht is equivalent to using a j^-difference which in 
turn is equivalent to the usage of a log-likelihood ratio. The 
numerator of this ratio is given by the likelihood ±fmodei of 
the model under study, e.g. the cMSSM. The denominator 
is given by the maximum of a phenomenological likelihood 
±fpheno which depends directly on the model predictions a,. 
These possess an expression as functions a, (P) of the model 
parameters P of -Z’modei- Note that in -Z^heno however, the 
a,- are treated as n independent parameters. We now identify 
- -21n±fmodei and x^^^ta - -21n±fpheno- When insert¬ 
ing a, (P), one is guaranteed to find 


.Z’pheno(ai (P),..., (P)) = ±fmodel (P) • 

Using these symbols, Zcombined written as 

^model(P) 


Zcombined 2 In 


°Zpheno(^l; ■ ■ ■ 


( 12 ) 


(13) 
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X^/ndf 


Fig. 15 Numerical example showing the distribution of jnAi using 
combined and split measurements. Using split measurements smaller 
values of ;(;^/ndf are obtained. Because in this example all measure¬ 
ments are Gaussian, this is equivalent to larger p-values. We call the 
effect of obtaining larger p-values when using split measurements the 
dilution of the p-value. 


where ai,... maximise -S^^heno- Note that in this formu¬ 
lation the model predictions a, do not necessarily need to 
correspond directly to measurements used in the fit, as it is 
the case for our example. For instance the model predictions 
Qi might contain cross sections and branching ratios which 
are constrained by rate measurements. 

Using ndfsplit = N, ndfcombined = « and ndfdata =N -n 
Eq. (pTI implies 


^.split 

ndfsplit 


y2 y2 

^combined Adata 


ndfcombined ^ ^ + tl 


(14) 


The more uncombined measurements are used, the larger 
N — n gets and the less the p-value depends on the first term 
on the right hand side, which measures the agreement be¬ 
tween data and model. At the same time, the p-value de¬ 
pends more on the second term on the right hand which mea¬ 
sures the agreement within the data. Especially, for n fixed 
and N 

^split _ Xdata 
ndfsplit ndfdata' 


Since in the case of purely statistical fluctuations of the 
split measurements around the combined value the agree¬ 
ment within the data is unlikely to be poor, the expectation 
is 


~ 1 (16) 

ndfdata 


even if there was a deviation between the model predic¬ 
tions and the physical combined observables. So most of the 
time the p-value will get larger when uncombined measure¬ 
ments are used, hiding deviations between model and data. 
As a numerical example Eig. 15 shows toy distributions of 


^combined ^5El!L fop one observable (n=l), ten measure- 

^^Icombined ^‘^fsplit 

ments (N=10) and a 3(7 deviation between the true value and 
the model prediction. We call this effect dilution of the p- 
value. It explains the large p-value for the Large Observable 
Set by the overall good agreement between the uncombined 
measurements. 

On the other hand if there is some tension within the 
data, which might in this hypothetical example be caused 
purely by statistical or experimental effects, the “innocent” 
model is punished for these internal inconsistencies of the 
data. This is observed here for the Medium Observable Set 
and Small Observable Set. Hence, and in order to incorpo¬ 
rate our assumption that ATLAS and CMS measured the 
same Higgs boson, we produce our own combination of cor¬ 
responding ATLAS and CMS Higgs mass and rate measure¬ 
ments. We also assume that custodial symmetry is preserved 
but do not assume that /z —> 77 is connected to h ^ WW 
and h —?> ZZ as in the official ATLAS combination used 
in Small Observable Set. We call the resulting observable 
set Combined Observable Set. Note that for simplicity we 
also combine channels for which the cMSSM model predic¬ 
tions might differ due to different efficiencies for the differ¬ 
ent Higgs production channels. This could be improved in 
a more rigorous treatment. Lor instance the could be de¬ 
fined by Eq. ( [Ts] ) using a likelihood .ifpheno which contains 
both the different Higgs production cross sections and the 
different Higgs branching ratios as free parameters a,. 

Using the Combined Observable Set we get a p-value of 
(8.3 ± 0.8)%, which is significantly smaller than the diluted 
p-value of (41.6 ±4.4)% for the Large Observable Set. The 
good agreement within the data now shows up in a small 
X^/ndf of 68.1/65 for the combination but no longer affects 
the p-value of the model fit. On the other hand the p-value 
for the Combined Observable Set is larger than the one for 
the Medium Observable Set, because the tension between 
the ATLAS and CMS measurements is not included. This 
tension can be quantified by producing an equivalent AT¬ 
LAS and CMS combination not from the Large Observable 
Set but from the Medium Observable Set giving a relatively 
bad x^/ndf of 10.9/6. 

Linally, as discussed briefly in Section we employ 
the Medium Observable Set again to compare the preferred 
parts of the parameter space according to the profile like¬ 
lihood technique (Lig. 2(b)| l with the parameter ranges that 
are preferred according to our pseudo fits. In Lig. |16(a) and 
|16(b)| we show the 1-dimensional distributions of the pseudo 
best fit values for Mq and Mi/2- The 68% and 95% CL re¬ 
gions according to the total pseudo best fit x^ are shown. 
As expected by the non-Gaussian behaviour of our fit, some 
differences between the results obtained by the profile like¬ 
lihood technique and the pseudo fit results can be observed. 
Lor the pseudo fits, in both parameters Mq and Mi)2, the 
95% CL range is slightly smaller compared to the allowed 
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Fig. 16 Distribution of the pseudo best fit values for 
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range according to the profile likelihood. Considering the 
fits of the pseudo datasets, Mq is limited to values <8.5 TeV 
and Mj /2 is limited to values <2.7 TeV, while the profile 
likelihood technique yields upper limits of 10 TeV and 3.5 
TeV, roughly. The differences are relatively small compared 
to the size of the preferred parameter space, and may well 
be an effect of the limited number of pseudo datasets that 
have been considered; the use of the profile likelihood tech¬ 
nique for the derivation of the preferred parameter space can 
therefore be considered to give reliable results. However, as 
discussed above, in order to get an accurate estimate for the 
95%-CL regions, the p-value would have to be evaluated at 
every single point in the parameter space. 


5 Conclusions 

In this paper we present what we consider the final analysis 
of the cMSSM in light of the LHC Run 1 with the program 
Fitting. 

In previous iterations of such a global analysis of the 
cMSSM, or any other more general SUSY model, the fo¬ 
cus was set on finding regions in parameter space that glob¬ 
ally agree with a certain set of measurements, either using 
frequentist or bayesian statistics. However, these analyses 
show that a constrained model such as the cMSSM has be¬ 
come rather trivial: because of the decoupling behaviour at 
sufficiently high SUSY mass scales the phenomenology re¬ 
sembles that of the SM with dark matter. This does how¬ 
ever not disqualify the cMSSM as a valid model of Nature. 
In addition, there are undeniable fine-tuning challenges, but 
also these do not statistically disqualify the model. There¬ 
fore, before abandoning the cMSSM, we apply one crucial 
test, which it has not been performed before: we calculate 
the p-value of the cMSSM through toy tests. 


A likelihood ratio test of the cMSSM against the SM 
would be meaningless, since the SM cannot acommodate 
dark matter. Thus we apply a goodness-of-fit test of the 
cMSSM. As in every likelihood test (also in likelihood ra¬ 
tio tests), the sensitivity of the test towards the validity of 
the model depends on the number of degrees of freedom 
in the observable set, while the sensitivity towards the pre¬ 
ferred parameter range does not. Thus, when calculating the 
p-value of the cMSSM, it is important that the observable 
set is chosen carefully. First, only such observables should 
be considered for which the cMSSM predictions are, in prin¬ 
ciple, sensitive to the choice of the model parameters, inde¬ 
pendent of the actually measured values of the observables. 
This excludes e.g. many LEP/SLD precision observables, 
for which the cMSSM by construction always predicts the 
SM value for any parameter value. Second, it is important 
that observables are combined whenever the corresponding 
cMSSM predictions are not independent. Otherwise the re¬ 
sulting p-value would be too large by construction. It should 
be noted that the allowed parameter space for all observable 
sets studied here is virtually identical. It is only the impact 
on the p-value which varies. 

In order to study this dependence, several observable 
sets are studied. The main challenge arises from the Higgs 
rate measurements. Since the cMSSM Higgs rate predic¬ 
tions are, in principle, very sensitive to the choice of model 
parameters, the corresponding measurements have to be in¬ 
cluded in a global fit. Using the preferred observable sets 
“combined” and “medium” (as described in Section [43] and 
Tab. [^, we calculate a p-value of the cMSSM of 4.9% 
and 8.4%, respectively. In addition, the cMSSM is excluded 
at the 98.7% CL if Higgs rate measurements are omitted. 
The main reason for these low p-values is the tension be¬ 
tween the direct sparticle search limits from the LHC and 
the measured value of the muon anomalous magnetic mo- 
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ment (g — 2)^. If e.g. (g — 2)^ is removed from the fit, the 
/7-value of the cMSSM increases to about 50%. However, 
there is no justification for arbitrarily removing one variable 
a posteriori. On the other hand, the observable sets could 
be artificially chosen to be too detailed, such that there are 
many measurements for which the model predictions cannot 
be varied individually. This is the case for the Large Observ¬ 
able Set of Higgs rate observables in Tab.[^ the inclusion of 
which does thus not represent a methodologically stringent 
test of the /7-value of the cMSSM. 

Thus, the main result of this analysis is that the cMSSM 
is excluded at least at the 90% CL for reasonable choices of 
the observable set. 

The best-fit point is in the f-coannihilation region at 
Mo « 500 GeV, with a secondary minimum in the focus- 
point region at Mi/ 2 ,Mo ^ 2TeV. A comparison of the p- 
values of coannihilation and focus-point regions can serve 
as an estimate of a likelihood-ratio test between a cMSSM 
at Mq « 500 GeV which can be tested at the LHC, and a 
“SM with dark matter” with squark and gluino masses be¬ 
yond about 5 TeV. Since the focus point manifests a more 
linear relation between observables and input parameters in 
the toy fits, and thus a more ^^-distribution like behaviour, 
it reaches a slightly higher /7- value than the f-coannihilation 
region. This shows that even the best-fit region offers no sta¬ 
tistically relevant advantage over the “SM with dark matter”. 
Thus, we can conclude that the cMSSM is not only excluded 
at the 90% to 95% CL, but that it is also statistically mostly 
indistinguishable from a hypothetical SM with dark matter. 

In addition to this main result, we apply the first com¬ 
plete scan of the possibility of the existence of charge 
or colour-breaking minima within a global fit of the 
cMSSM. In addition, we calculate the lifetime of the best 
fit points. We find that the focus-point best-fit-region is sta¬ 
ble, while the f-coannihilation best-fit region is either stable 
or metastable, with a lifetime significantly longer than the 
age of the Universe. 

It is important to note that the exclusion of the cMSSM 
at the 90% CL or more does in general not apply to less re¬ 
stricted SUSY models. The combination of measurements 
requiring low slepton and gaugino mass scales, such as 
(g — 2)^, and the high mass scales preferred by the SM- 
like Higgs and the non-observation of coloured sparticles 
at the LHC puts the cMSSM under extreme tension. In the 
cMSSM these mass scales are connected. A more general 
SUSY model where these scales are decoupled, and prefer¬ 
ably also with a complete decoupling of the third generation 
sleptons and squarks from the first and second generation, 
would easily circumvent this tension. 

Therefore, the future of SUSY searches at the LHC 
should emphasize the coverage of any phenomenological 
scenario which allows sleptons, and preferably also third 
generation squarks, to remain light, while the other sparti¬ 


cles can become heavy. Many loopholes with light SUSY 
states still exist, as analyses as in | |159| show, and there exist 
potentially promising experimental anomalies which could 
be explained by more general SUSY models [160 - 1621. 

On the other hand, the analysis presented here shows that 
SUSY does not directly point towards a non-SM-like light 
Higgs boson. The uncertainty on the predictions of ratios of 
partial decay widths and other observables at the LHC are 
significantly smaller than the direct uncertainty of the LHC 
Higgs rate measurements. This is because of the high SUSY 
mass scale, also for third generation squarks, imposed by 
the combination of the cMSSM and the direct SUSY parti¬ 
cle search limits. These do not allow the model to vary the 
light Higgs boson properties sufficiently to make use of the 
experimental uncertainty in the Higgs rate measurements. 
This might change for a more general SUSY model, but 
there is no direct hint in this direction. The predicted level 
of deviation of the light Higgs boson properties from the 
SM prediction at the ^(1%) level is not accessible even at a 
high-luminosity LHC and requires an collider. 

In summary, we find that the undeniable freedom in 
choosing the observable set - before looking at the exper¬ 
imental values of the results - introduces a remaining soft¬ 
ness into the exclusion of the cMSSM. Therefore, while we 
might have preferred to find SUSY experimentally, we find 
that at least we can almost complete the second most revered 
task of a physics measurement; with the combination of as- 
trophysical, precision collider and energy frontier measure¬ 
ments in a global frequentist analysis we (softly) kill the 
cMSSM. 
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